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(57) ABSTRACT 

Linear space-time multiuser detectors are described that 
combine array processing, rake detection, and multiuser 
detection and that satisfy a minimum mean-squared error 
criteria. Both embodiments can be implemented adaptively 
so as to account for unknown sources of interference. The 
adaptive implementation of the first embodiment of a linear 
space-time detector requires explicit estimates of the chan- 
nel coefficients. On the other hand, the adaptive implemen- 
tation of the second embodiment does not require these 
estimates since they can be adaptively obtained. Both of 
these detectors provide significant performance advantages 
over single-user space-time detectors and over multiuser 
detectors which do not account for CDMA interference from 
unknown sources. 
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DETECTORS FOR CDMA SYSTEMS 

FIELD OF THE INVENTION 

This invention relates generally to direct sequence (DS) 
code division multiplex access (CDMA) systems subjected 
to multipath fading and multiple access interference (MAI), 
and more particularly to signal detection techniques for use 
in such systems. 

BACKGROUND OF THE INVENTION 

In wireless CDMA systems such as those proposed for 
next generation mobile wireless standards, the major 
impediments at the physical layer are multipath fading and 
multiple access interference caused by co-channel users 
which are not orthogonal to the desired user. The rake 
receiver, described in R. Price and R E. Green, "A Com- 
munication Technique for Multipath Channels," Proceed- 
ings of the IRE, Vol. 46, pp. 555-570, March 1958, attempts 
to combat multipath fading by coherently combining resolv- 
able multipath replicas of the desired signal. Multiuser 
detection described in S. Verdu, "Multiuser Detection/' 
Cambridge University Press, New York, 1998, addresses the 
problem of MAI by actively accounting for its presence 
when detecting the desired user. 

More recently, there has been growing interest in using 
array processing for further improving receiver perfor- 
mance. These techniques have focused on using multiple 
antennas at the base station receiver to provide antenna gain 
and/or diversity gain and allow the possibility of spatial 
processing. By combining these space-domain techniques 
with time domain techniques like rake detection and mul- 
tiuser detection, the resulting space-time detectors show 
promise of improving the capacity of CDMA systems as 
compared to traditional time-domain-only detectors. See, for 
example, A. Paulraj and C. Papadias, "Space-Time Process- 
ing for Wireless Communications," IEEE Signal Processing 
Magazine, Vol. 14, No. 6, pp. 49-83, November 1997. The 
first generation of space-time CDMA detectors used array 
processing with either rake detection or multi-user detection. 
See, respectively, A. Naguib and A. Paulraj, "Performance of 
Wireless CDMA with M-ary Orthogonal Modulation and 
Cell Site Antenna Arrays," IEEE Journal on Selected Areas 
in Communications, Vol. 14, No. 9, pp. 1770-1783, Decem- 
ber 1996 or S. Miller and S. Schwartz, "Integrated Spatial- 
Temporal Detectors for Asynchronous Gaussian Multiple - 
Access Channels," IEEE Transactions on Communications, 
Vol. 43, No. 2/3/4, pp. 396-411, February/March/April 
1995. Later space-time CDMA detectors combined all three 
processing techniques. See: H. Huang; S. Schwartz, S. Verd 
u, "Combined Multipath and Spatial Resolution for Mul- 
tiuser Detection: Potentials and Problems," Proceedings of 
the IEEE International Symposium on Information Theory, 
p. 380, 1995; or M. Nagatsuka and R. Kohno, "A Spatially 
and Temporally Optimal Multi-User Receiver Using an 
Array Antenna for DS/CDMA," IEICE Transactions on 
Communications, Vol. E78-B, No. 11, pp. 1489-1497, 
November 1995. 

While the foregoing systems operate satisfactorily, 
improvements can be made, in particular to the space-time 
detectors which combine all three processing techniques. 
While the detector in the Nagatsuka and Kohno paper is 
optimum in the maximum likelihood sense, its computa- 
tional complexity is exponential with respect to the number 
of users. Hence it is too complex to implement for practical 
systems. In the paper by Huang, Schwartz and Verdu, a 
tradeoff between performance and complexity is made, but 
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this detector was not implemented adaptively since it used a 
zero-forcing criteria. Adaptive implementations allow 
receivers to account for unknown sources of interference 
thus improving the detector performance and increasing the 
5 system capacity. For example, a base station receiver could 
account for interference from adjacent cells or from an 
embedded microcell, while a handset receiver could account 
for interference from signals it is not explicitly demodulat- 
ing. 

10 SUMMARY OF THE INVENTION 

This invention detects DS-CDMA signals utilizing a rake 
receiver, array processing, and multiuser detection. When 
combined with array processing, the rake receiver is often 

15 called a space -time rake receiver. The invention uses a 
minimum mean-squared error (MMSE) criterion in the 
multiuser detector. This criterion allows for relatively simple 
implementation in the form of a linear detector and also 
allows for adaptive implementation. Adaptive implementa- 

20 tions are useful in practical situations where there is limited 
knowledge of the various received signals. As discussed 
above, both uplink and downlink capacity can be improved 
using adaptive detectors which account for unknown inter- 
ference. 

25 For pulse amplitude modulated (PAM) data signals, two 
embodiments of the invention provide options for trading off 
between performance and adaptive implementation com- 
plexity. The first embodiment performs better if there is 
perfect knowledge of the users' signal parameters. On the 

30 other hand, it requires more explicit channel information for 
adaptive implementation. Suppose the receiver consists of 
P^l antennas and the received signal consists of K 
DS-CDMA data signals, each with L delayed/weighted 
multipath replicas. At each antenna, a bank of filters is 

35 matched to the KL spreading codes with their multipath 
timing delays. The filter outputs are weighted according to 
the complex conjugate of the estimated channel (multipath 
and array) parameters and combined to form a K-vector 
where each component corresponds to one of the K codes. 

40 The real part of each component is taken. Estimates of the 
multipath delays and channel parameters may be obtained, 
for example, from a training or pilot signal. The conven- 
tional space-time rake receiver would pass each vector 
component to a decision device for estimating the corre- 

45 sponding user's PAM data symbol. However, unless the 
users* signals are orthogonal in the space-time domain, the 
vector components are contaminated with multiaccess inter- 
ference from the other users. This invention proposes the use 
of a linear combiner prior to the decision devices to suppress 

50 the multiaccess interference. This linear combiner is repre- 
sented by a real K-by-K matrix W A which multiplies the real 
K-vector. The matrix minimizes the mean squared error 
between its product with the real K-vector and the K-vector 
of data symbols. Because the matrix uses the minimum 

55 mean -squared error (MMSE) criterion, well-known adaptive 
algorithms can be used to adaptively obtain it. Each com- 
ponent of the final K-vector output corresponds to one of the 
K codes and is passed to either a decision slicer or a decoder 
for further processing. 

60 The second embodiment has marginally inferior perfor- 
mance when the channel parameters are known exactly. 
However, under practical conditions of channel mismatch, 
the second embodiment often performs better. In terms of 
adaptive implementation, this embodiment requires less 

65 information (it does not require explicit channel estimates) 
but may be slower to adapt. As with the first embodiment, 
the front end of the second embodiment consists again of a 
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bank of KL matched filters for each of the P antennas sequences s^t). For Detector A, the data sequence is pulse 
followed by weighting and combining. However, unlike the amplitude modulated (PAM); for Detector B, the data 
first embodiment, the order of the real operator and linear sequence can be either PAM or quadrature PAM (QAM), 
combiner (K-by-K matrix multiplication) are exchanged. In Jhe transmitted signal undergoes frequency-selective fading 
doing so, the weighting, combining, and linear combiner can s in the channel and arrives at the receiver as L time- 
be represented by a single complex K-by-KLP matrix W^. resolvable multipath components with complex fading chan- 
The real parts of the components in the resulting K-vector nel coefficients c^t) . . . c kfL (t). Assume that the receiver is 
arc passed to either a decision device or a decoder for further a p - e l ement lincar arrav - If each resolvable multipath corn- 
processing. In the adaptive implementation, adaptive algo- P oaent a P lanar ^vefront with angle G^t) wUh 
•*l u j * u* • *u * * wt l. u i „ respect to a linear array and if the array spacing is sum- 
nthms can be used to obtain the matrix W»; hence channel 10 . , . , . « / t t , • j? . i . ■ 

A . j ciently close (e.g., A./2 ) so that there is perfect correlation 

estimates are not explicitly required. ^ % afray £ f & giyen wavefr £ nt) ^ ^ phase 

For quadrature PAM (QAM) data signals, real operators oflfeet of me pth element, with respect to the first is h kl 

are not required in the detector since the signal constellation (t)=exp[[n(p-l) sin B k fi)]]. Trie received signal at the pui 

is two-dimensional. In this case, a variation on the second antenna for given symbol period (and ignoring intersymbol 

embodiment would be used which does not use the real 15 interference) is: 
operators but follows the linear combiner directly with the 

appropriate decision device or decoder. K L (ij 

While the MMSE detection techniques are powerful in M') = £ Z W0c w M*(r-r u M*fr('-r w )+n p (O 

and of themselves, their performance can be further t=1 ' =1 

enhanced by using them with other multiuser detection 20 

techniques such as interference cancellation where where A^ is the amplitude for user k, x kl is the delay for the 

interference, instead of being projected away as is done with fcth user's 1th multipath, and n^t) is the additive Gaussian 

linear multiuser detectors, is explicitly subtracted from the noise process which accounts for out-of-cell interference 

received signal. Interference cancellation can occur prior to and background noise. We now make the following assump- 

or after the MMSE linear combiner. 25 tions to simplify the analysis. 

Both embodiments of the MMSE detector can be gener- 1. The signals are received bit synchronously, 

alized for practical systems which utilize, in addition to the 2. The time spread is small compared to the symbol period 

data bearing channels, supplemental channels which may so that intersymbol interference can be ignored, 

serve as pilot signals. Also, the MMSE detectors can operate 3. The phase offsets and channel coefficients are constant 

in systems where signals are transmitted with different 30 over a symbol period. 

spreading factors. The first assumption will be dropped later. Under these 

A CDMA signal detector in accordance with the invention assumptions, the chip matched-filter output for the received 

can provide significantly improved performance relative to si 8 nal m equation (1) is a complex N-vector 

a conventional signal detector which does not use multiuser 35 

detection or which does use multiuser detection but does not A 

adaptively mitigate unknown interference. tp = £ L Wm*u^ +», 

BRIEF DESCRIPTION OF THE DRAWINGS 

The foregoing and other features of the invention may wne re s kt is the chip matched-filter N-vector corresponding 

become more apparent when the ensuing description is read 40 t 0 s^t-x* ,) and n p is the complex N-vector corresponding to 

together with the drawings, in which: the Gaussian noise. The spreading codes are normalized to 

FIG. 1 shows a prior art single-user space-time rake have unit energy: ||s^||-l. The spreading codes are assumed 

receiver employing a multi-element antenna array, a bank of to be random. However, adaptive implementations of these 

N-chip filters at each antenna, and a channel weighter and ^ MMSE detectors require the use of short spreading codes 

combiner; which repeat after a few symbol periods. The following 

FIG. 2 shows the generic embodiment of the invention; notation will be used hereinafter: 

FIG. 3 shows a first embodiment (Detector A) of the 

multiuser space-time MMSE detector; 

FIG. 4 shows a second embodiment (Detector B) of the r n u n. i. tr i « > r 

, . . * h k i «• [hw. . . hkj. p] complex P-vector of array 

multiuser space-time MMSE detector; * coefficients 

FIG. 5 shows an adaptive implementation of the second H - [h^ . . . \^ . . . h V(1 . . . hjcd complex P x kl array matrix 

embodiment; RD = ^SCV • • ■ H JtL . . . . . . complex KLP x KL array matrix 

FIG. 6 shows interference cancellation located prior to the ■=■ [c^ . . . c^ J 1 * complex L- vector of channel 

linear combiner of the first embodiment; 55 coeffic'ts 

FIG. 7 shows interference cancellation located prior to toe £ - £ * x E^Sf 

linear combiner of the second embodiment; and matrix 

FIG 8 shows interference cancellation following the defined by R^-o^^-OK+ii) 

linear combiner of the multiuser space-time MMSE detector J . , . " s 1 ^ Sk " J ,. J 

, . . . . 4 . „ A = diagfA. . . . Ajj real K x K amplitude matrix 

in accordance with the invention. 60 b - [b, b K f real K-vector of data bit. 

FIG. 9 shows a block diagram of a post-combiner inter- K uxu identity matrix 

ference cancellation receiver in accordance with the inven- \» u-vector of ones 

tion. 



u x u matrix of zeros 
0 U u-vector of zeros 



GENERAL DESCRIPTION 



65 

Consider a K user system where the kth user (k=l . . . K) The corresponding estimated values for the array and chan- 
modulates its data sequence bj^t) with an N-chip spreading nel coefficients will be represented with the symbol over 



07/13/2004, EAST Version: 1.4.1 



receiver is: 



US 6,301,293 Bl 
5 6 

the symbol for the value. The noise vector is a zero-mean 

complex Gaussian vector whose distribution can be written IRetf)] fro*! f Re{R) /m(ft)T) 

in terms of its (component-wise) real' and imaginary com- tm(ff) ~i } <^ 

ponents: v L K * y )V 

5 

rft*01 rro^l J Im Using ^ &ct that (ft^RH^C-C'Kft^oRlC 

[/m(n)J VL°Af J [Qnxn In \) (where o denotes the component- wise product between two 

same -sized matrices), the K vector at the output of the 
where we define the real and imaginary operations for 10 com ^ers and the input to the decision devices 

matrices and vectors as Re(X)=(X+X*)/2 and Im(X)=(X- 14 *" 14 * can be concisely written as: 
X*)/2 where * denotes the complex conjugate. Hence Re(n) ^ j-R e {(ft D C:) w z }MAb+n c , jrfV (3) 

and Im(n) are zero-mean Gaussian random vectors whose ccnv 

components have variance o 2 and are mutually uncorrelated. where H denotes the complex conjugate transpose, o rep- 

^ . _ _ , v 15 resents the component-wise product of two same-sized 

Conventional (Space-Time Rake) matrices, and M=Re{e"[(ft*)oR]C}. The noise vector is 

Receiver, (FIG. 1) purdy reaJ and hag ^^0^ 

In the context of this array-multipath channel, the prior art . 
single user detector is a correlator matched to the desired ^^(^o™)- 

user's composite array-multipatb-spreading code signal. 20 ^ fy^f^ffy^ If thc data fa BPSK modulated> 

This detector does not account for the presence of interfer- , he w , fof ^ > k is si , ^ hanJ ^ of ^ k(h 

ers; however, it is the maximum-likelihood detector if there component of Rely V b =sen?Re<v V) The corre 
are no interferes or if they are orthogonal to the desired user dj bft ^ p 0 > the "conventional rake 

m array-code space. As shown in FIG. 1, the detector 
consists of a bank of correlators 10 at each antenna, matched 25 
to the KL multipath spreading codes s hl . . . s^ . . . s K1 . . . 
s KiL . The notation <s w ,-> in the left- most boxes of FIG. 1 p^V) - Q 

indicates taking the inner product between s M (the kth user's * 
spreading code corresponding to the 1th multipath delay) and 
the input vector to the box. The dot in the notation <s k y> 30 

represents the input. The timing estimates are obtained using where is the (k,k)th element of matrix X. The perfor- 

a separate timing estimation algorithm. The output of the mance comparisons among the various detectors are given in 

upper left-hand box (which corresponds to the first compo- terms of BPSK modulated data, but they can also be done for 

nent of the vector z 7 ) is <s x 1 jff ,r 1 >=s 1 a h t 1 where the H any QAM modulated data in general, 
superscript signifies the Hennitian transpose (take the com- 35 

plex conjugate of each element, then take the transpose of Generic Space-Time Linear Multiuser Rake 

the resulting vector) of a complex vector. The matched filter Receiver (FIG. 2) 

outputs for at least a subset of the multiple signals, a ™ . ^ , • . j t_ • 

i ... c u . . r • i j The prior art space -tune rake receiver descnbed above is 

plurality of the multipath components of these signals, and . * ■ • ,i_ «t_ . c • . L 

! 1t , r . ■ j- A j ■ *t_ j • An a single user receiver in the sense that for a given user, the 

a plurality of receiver antennas are indicated in the drawing 4 ° , , . c , f * o- ■* 

at (z z ) demodulation uses information only from that user. Since it 

^ 1 ' ' ' p ' does not account for the presence of interference from other 

At reference number 11, the correlator outputs are uscrSj its performance suffers . FIG. 2 shows a generic 

weighted by the complex conjugate of their corresponding mu i thlser spa ce-time rake receiver. It accounts for the mul- 

channel estimate. More specifically, the correlator output at Uaccess interference usmg a combiner. Two versions 

the pth antenna for the 1th multipath of user k is weighted by of ^ detector are giverK first> ^ detector A, 

the complex conjugate of the estimate of the corresponding u$es the real operators before me linear com biner. In this 

channel (array/multipath) coefficient h WjP c^ kV These esti- casej re f erence 2 0 in FIG. 2 consists of a bank of real 

mates are obtained using a separate channel estimation oper ators followed by a linear combiner, and each of the 

algorithm. At each antenna, the L components for the kth componerjts i4l . . . 14* is a sheer which determines the 

user are added 12), and then the resulting P components for closest estimated symbol to its input. The second, detector 

the kth user are added 13* (k-1 . . . K). Each component is B> ^ a combincr represented by a complex matrix 

then passed to either a decision device or a decoder. For followe d by the real operators. In this case, reference 20 

PAM data signals, the decision device for the kth user 14„ COQsists only of ^ combincr) arjd ^ decision 

k-1 . . . K) gives as its output the symbol which is closest devices consists of a real operator and sUccn Despite lhe 

in the Euclidean distance sense to the realpart of akthvector simi iarities, the fact that the linear combiners for detectors 

component. For QAM data signals, the real operator is not A and B results ^ different performances and adaptive 

required. implementations. We now describe the two detector embodi- 

From FIG. 1, the KLP-vector z=[z/* . . . Z/] T at the men ts more thoroughly. 



output of the matched filters can be written as: 

z=RH D CAb+n* (2) 



60 



Space-Time Linear Multiuser Detector A (FIG. 3) 



A set of sufficient statistics for the K users is given in 
where _ we define the KLPxKLP block-toeplitz equation (3) as the K-vector Re{y ce>m ,}. The objective of the 
matrix RoRQI^p and tP represents the KLPxKL matrix, linear MMSE detector is to apply a KxK linear combiner to 
Q denotes the Kronecker product operation between two 65 this vector such that the mean-squared error between its 
matrices, and n s denotes the complex noise vector is a resulting vector and the data vector b is minimized. In other 
complex Gaussian random vector with distribution words, find the real KxK matrix W A such that 
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W A =argmin{^Q|WRc{y £ « IV }-frt| 2 ]}. 



W^-Etb(Re{y^"}) w ]{E[Re{y« :o " t '}Re w {y w ' 1 '}]}- 1 
^AM(MA 2 M"+o 2 lk)- 1 
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definite and hence invertible; and it follows that RejC" 
(5) [(H^oRJC} is also invertible. 

Let us now calculate the bit error rate of this first MMSE 
detector (MMSE A). The K-vector at the slicer inputs is: 



This detector in essence strikes the proper balance in com- 
bating both residual Gaussian noise, manifested through 
n con»> a nd multiple access interference, manifested through 
the off-diagonal terms of Re{C'T;(ft /f H)oR]C}. It can be 
shown that the solution is: 



(6) 



FIG. 3 shows a block diagram of detector-Ain accordance 
with the invention. It contains P antennas and a bank of 
matched filters 10. The box 31 represents a linear combiner 
defined by the K-by-K matrix W A given by equation (6). The 
matrix can be put into words as follows: it is the product of 
3 terms: (a) the diagonal amplitude matrix (A); (b) the 
equivalent array/channel/code correlation matrix (M); and 
the matrix inverse of the matrix (MA^^+^M) which 
consists of the product of the array/channel/code correlation 
matrix (M), the square of the amplitude matrix (A 2 ), and the 
Hermitian conjugate of the array/channel/code correlation 
matrix (M 7 *) summed with the estimated array/channel/code 
correlation matrix (M) which is weighted by the background 
noise variance (cr 2 ). This matrix minimizes the mean- 
squared error criterion given in equation (5), and because of 
the convexity of the mean-squared error expression in (5), 
the MMSE solution (6) represents a global minimum. A 
bank of decision slicers (32^32^ follows the linear com- 
biner. 

An attractive feature of the MMSE detector is that it can 
be implemented adaptively using well-known adaptive algo- 
rithms (33) such as least-mean-squares or recursive least- 
squares. Using a training signal (for example, the user data 
bits), these techniques can be used to adaptively obtain W A . 
Adaptive implementation is an option if direct calculation of 
W A is deemed too complex. Note that the adaptive imple- 
mentation requires knowledge of the K users* spreading 
codes, multipath delays, and channel parameters. 

Assuming that the array and channel estimates are exact 
C=C; ft=H, we can rewrite the expression for W A in (6) as: 



W^-A-^Rc^fH^oRlCj+c^A- 2 ]- 1 , 

and another attractive feature is that MMSE detector 
approaches the decorrelating (zero-forcing) detector (up to a 50 
scalar factor) as the residual Gaussian noise approaches 
zero: 



In other words, from (3) we have that: 



35 



Km{W A ) « A- L (Re[C H [{H H H)cR]C})' 1 . C) 



where 



10 



Hence the bit error rate for the kth user using the MMSE-A 
detector is: 



15 



20 



25 



(W A MAb) k 



(8) 



Letting W^A-^RefC^^HJoRjC})- 1 and assuming per- 
fect channel estimates, the bit error rate for DD-A, is: 



r^[{Re{C"[(H»H)°R]C)r% M 



(9) 



One drawback of this MMSE detector is that the estima- 
tion of the channel and array coefficients cannot be incor- 
porated into the adaptive algorithm for obtaining W A . The 
estimates must be explicitly obtained using some other 
means (like a training or pilot signal). This leads us to ask, 
"is it possible to design an adaptive MMSE which does not 
require prior knowledge of the array and channel param- 
eters?" The answer is a "yes," as we will see in the following 
section. 

Space-Time Linear Multiuser Detector B (FIG. 4) 

The goal of this subsection is to derive an MMSE detector 
which can be implemented adaptively and which does not 
require explicit array and channel estimates. (It will, 
however, still require timing estimates for the multipath 
delays.) FIG. 4 shows the basic structure of Detector B 
which exchanges the places of the real operators and linear 
combiner of Detector A. Following this exchange, the result- 
ing cascade of the array combining, multipath combining, 
and linear combiner can be collected into a single complex 
KxKLP matrix multiplication called W fl , as shown in FIG. 
5. Specifically, we define this matrix using the MMSE 
criterion as: 

W B =zr & mn(E[\\Wz-ty\ 2 ]} t 



55 



lim(H' /1 y)=A- 1 (Re|C w [(H w //)o/fiq)" 1 /fc{C ,/ [(// w //)o/fjCM* 

(T-»0 g0 

= b 



and it follows from (4) that 

W B = E(b?)[E(zz H )T l 

= AG h R[RGA 2 G h /? + cr^ff] -1 
= AG H \RGA 2 G H + o^/iop]" 1 



(10) 



and the multiple access interference is forced to be zero. Id where GoH^C is the array/channel matrix. FIG. 5 shows a 

the implementation of the decorrelating detector A (DD-A), 65 block diagram of the MMSE-B detector in accordance with 

we replace W A given by equation (6) with W A aA^Re" the invention. It contains P antennas and a bank of matched 

i{M}. We will assume that C/tXH^HJoRjC is positive filters 10. The box 50 represents a linear combiner defined 
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by the K-by-KLP matrix given by equation (10). This 
matrix minimizes the mean-squared error criterion given 
above by 

W B =Brffnm\E[\\W Z - bf}}. 



The last line of equation (10) states that the matrix is the 
product of 3 terms: (a) the diagonal amplitude matrix (A); 
(b) the array/channel combiner matrix H D C which is equiva- 
lent to the array/channel combiner in the conventional 
detector but using the actual parameters and not their 
estimates; and (c) the matrix inverse of the matrix 
RGA 2 GH+o 2 I J<Li ». The matrix inverse consists of the prod- 
uct of the correlation matrix (R) defined following equation 
(2), the array/channel combiner matrix G, the square of the 
amplitude matrix (A 2 ), and the Hermitian conjugate of the 
array/channel combiner matrix summed with the KLP-by- 
KLP identity matrix weighted by the background noise 
variance (cr 2 ). If the data is PAM., a bank of real operators 
(3(^-30^) and decision slicers (32^32^) follows the linear 
combiner. Otherwise, if the data is QAM, only the complex 
decision slicers (1^-14^) are required. For PAM data, 
Detectors A and B, while similar in structure, are not 
equivalent due to the different placement of the real operator. 

To calculate the bit error rate of Detector B for BPSK 
data, the K-vector at the slicer inputs is 



where 



[Re{W B Rti D c)Ab) k 



By concatenating the P received N-vectors into a single 
NP-vector x^T 1 r r 2 r . . . r^ 7 ] 7 , we can combine the opera- 
tions of the matched filter bank with the linear combiner W B 
to create a single K-by-NP matrix which multiplies r. As 
shown in FIG. 6, the detector structure for the kth user can 
be reduced to the kth row of this K-by-NP matrix 60 
followed by the real operator 62 and decision slicer 63. 
Hence the MMSE-B detector can be implemented with a 
simple tap-weight filter architecture. Its corresponding deco- 
rrelating detector (described below), the MMSE-A detector 
and the DD-A detector can likewise be implemented with 
this architecture. Like the MMSE-A detector, this version of 
the MMSE-B detector can be implemented adaptively using 
well-known adaptive algorithms (61). However, MMSE-B 
has the advantage of not requiring the array and channel 
estimates. The tradeoff is that the adaptation may be slower 
because there are more taps to adjust. 

While direct analytic comparisons between the two detec- 
tors is difficult, we will show that under perfect channel and 
array estimates, the performance of MMSE-A is uniformly 
superior to that of MMSE-B as the noise floor approaches 
zero. First, we show that the corresponding decorrelating 
detector for MMSE-B, DD-B, simply replaces W B with 



10 



Hence the bit error rate for the kth user using the MMSE-B 
detector is: 



(11) 



50 
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A-^fH^q^RH^CJ-^^q". (Recall that we have already 
made the assumption that (H z> C/ f RH i> C-C ff [(H^H)oR]C is 
invertible.) 

PROPOSITION 1 

hm[{A(H D C)f[RH D (CA 2 {H D C)f +o 2 I KLP \' 1 ] « 

A- l [{H D Cf%H D C\-\H D C) H 

Proof: Multiplying both sides by RH^CA^q^+o 2 !^ 
taking the limit as a approaches 0, we have that 

35 A(H D C) H =^A- l \(H D Q H RH D C]~\H D Q H [RH D CA 2 (H D cf + 



: A(H D C) H + \im(a- 2 A- l [(H D Q H RH D C]~ l (H D C) H ) 



20 



Hence, from the fact that (H^q^RH^C-C^^oRjC 
25 and the fact that the array and multipath combiners in 
MMSE-A can be represented by (H D D) M , the DD-B is the 
bank of matched filters, followed by the combiners, fol- 
lowed by the decorrelator (C^ftH^HJoRlC)" 1 , followed by 
the real operators, and followed by the slicers. The DD-A is 
3 q nearly identical except that the order of the decorrelator and 
the real operators are switched, and the decorrelator is 
(Re{C H [(H /f H)oR]C})- 1 . By replacing with U^A' 1 
[(J^CfifoPCY^Cf in (11), the bit error rate of DD-B 
is: 



35 



<ry]Re[{C»[{H»H)oR]C)-% 



(12) 



We would now like to show that DD-A has a uniformly 
lower bit error rate than DD-B. 

PROPOSITION 2 



Define MeC^t^HJoRJC, which is Hermitian and 
assumed to be positive definite. It follows that [Re(M)]"\ is 
also positive definite, and it is sufficient to show that the kth 
diagonal element of [Re(M)]~ 1 is greater than or equal to the 
kth diagonal element of M~\ M~\ (Note that since M is 
Hermitian, its inverse is also Hermitian, hence the diagonal 
elements of the inverse are real.) Chose y to be real and set 
x4Re(M)] _1 y, and note that x is real. Thus x T1 ^ xT 
Re^x-yCRe^l^y and also x^y-y^Re^]" 1 ^ Using a 
generalization of Bergstrom's inequality which states that 
for any real x, real y, and positive definite M, x 7 MXxy r 
(M^y^x 7 }') 2 we have that 

and the result follows by using the kth unit vector for y. 
Despite this relationship, it turns out that the performance of 
DD-A and DD-B are actually quite similar. More 
65 importantly, it has been shown that the performance of 
MMSE-B can be better than MMSE-A for the practical case 
of inexact array and channel estimates. 



60 
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Extensions of MMSE Detectors to the . 
Asynchronous Multirate Case 

So far, we have assumed bit synchronous operation of 
these linear multiuser detectors. We now show how their 
operation can be extended to the bit asynchronous case. 
First, we consider the operation of a decorrelating detector 
in a synchronous single -rate environment. Similar reasoning 
applies to MMSE detectors; however, for pedagogical 
purposes, we focus on the decorrelating detector. Suppose 
that there are two users with respective spreading codes s 1 (t) 
and s^t). Table A below shows the relationship between the 
two codes in time, with the thin vertical lines representing 
symbol boundaries. The received signal is r(t)=A 1 s 1 (t)b,+ 
AaS^tJbj+nfl) where A* is the kth user's amplitude, b* is the 
kth user's data bit, and n(t) is the additive Gaussian noise. 

TABLE A 



2 user synchronous single-rate system 











■2(0 





TABLE B 



2 user asynchronous single-rate system 





Ol(0 






■3(0 


*2<t-T) • 



TABLE C 



2 user synchronous multi-rate system 











82(0 


sj(t-T) 


s 2 (t-2T) 





If user 1 is the desired user, its decorrelating detector for 
a given symbol interval is the matched filter for code Sj(t) 
projected into the null-space of s 2 (t). For the asynchronous 
single-rate case, shown in Table — B, the decorrelating 
detector for user 1 is projected into the null-space spanned 
by the linear combinations of user 2's codes which overlap 
during that symbol interval. If we let S2*R(t) be the portion 
of s^t) which overlaps with s^t) and let s^t-T) be the 
portion of s^t-T) which overlaps with s^t), men, assuming 
BPSK data modulation, the null space is spanned by s*(t)+ 
s/O-T) and s 2 ii (t)-s 2 z '(t-T). This idea can be extended to 
the asynchronous multi-rate case shown in Table — C. Using 
similar function definitions, the decorrelator for user 1 ties 
in the null space spanned by s^t^s^t-TJ+s^Ct-lT), 
(d+SaO-TJ-s^t-ZI), s 2 *(t)-s 2 (t-T>s/(t-2T), and s* 
(t)-s 2 (t-T)-s 2 / '(t-2T). As the disparity between data rates 
increases, the effectiveness of a linear multiuser detector 
decreases since it will be constrained to a more restrictive 
subspace. The reasoning above can be applied to MMSE 
detectors by similarly combining the appropriate subspaces 
of the interfering waveforms. 

MMSE Detectors and Interference Cancellation 
The MMSE techniques described here can be used in 
conjunction with a non-linear multiuser detection technique 
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known as interference cancellation. Whereas the linear mul- 
tiuser detectors rely on subspace projections to mitigate 
interference, the interference cancellers directly subtract off 
the interference. Two types of interference cancellation can 

S be used: pre-combiner, where the cancellation occurs before 
the linear combiner, and post-combiner, where it occurs 
after. Pre-combiner cancellation is used for canceling inter- 
ference which is not accounted for by the linear combiner. 
For example, we may choose to cancel a high-powered, 

10 high-data rate signal from within the cell instead of account- 
ing for it with the linear combiner. Post-combiner cancella- 
tion on the other hand is used for refining the symbol 
estimates made by the linear combiner. Using preliminary 
symbol estimates from the linear multiuser detector, the 

15 signals are reconstructed and subtracted from the received 
signals to form enhanced signals which are used for a second 
stage of symbol estimation. In both interference cancellation 
techniques, capacity can be potentially increased by sup- 
pressing the interference. 

20 FIG. 7 shows a block diagram of a pre-combiner inter- 
ference canceller in accordance with Detector A shown in 
FIG. 3. At the kth input to the linear combiner W^, there is 
interference due to users j«l . . . K, j*k, which will be 
accounted for with W A . However, there may also be other- 

25 wise unaccounted interference, for example, high-powered 
interference from within the cell. This interference could be 
subtracted from the input to if the spreading codes, 
multipath delays, data bits, and channel parameters at all 
antennas are known. In this case, the baseband signal for the 

30 interference could be reconstructed and its contribution to 
the kth input of W A could be estimated (by passing the 
reconstructed signal through the processing chain for user k) 
and subtracted at the point designated 70^=1 . . . K) 
Contributions from other interference for this or other inputs 

35 can be similarly computed and subtracted. 

FIG. 8 shows a block diagram of a pre-combiner inter- 
ference canceller in accordance with Detector B from FIG. 
4. Using the same method as described above for the 
Detector A pre-combiner interference canceller, the interfer- 

40 ence which is not accounted for by W B can be subtracted 
from each of its inputs at the point designated 80. 

FIG. 9 shows a block diagram of a post -combiner inter- 
ference cancellation receiver in accordance with the inven- 
tion. This receiver includes P antennas and has three stages, 
designated 90, 91 and 92. The first stage 90 receives 
baseband signals t p and makes preliminary symbol estimates 
for all K users with either a MMSE-A detector of FIG. 1 or 
a MMSE-B detector of FIG. 2. The second stage 91 uses the 
preliminary symbol estimates from the first stage 90 and 
knowledge of the K users' spreading codes, delays, and 
channel parameters to reconstruct the baseband received 
signals for each of the users. The second stage 91 subtracts 
the multiaccess interference with respect to a desired user k 

55 from the received signals r 3 . . . r P to form the following 
enhanced received signals for user k: 

K L 

60 

Note that the multipath interference with respect to user k 
can also be removed. We will not describe this option in 
detail because its relative gains will generally be negligible 
65 unless the number of users K is very small. Ideally, if the 
preliminary symbol estimates and the channel estimates are 
perfect, there will be no multiaccess interference in the 
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enhanced received signals. However, this will generally not means for processing the output of the linear combiner 
be the case in practice. In the third stage 82, the enhance to produce estimates of the respective data symbols, 
received signals are processed with a conventional space- 5. A method of demodulating K direct sequence code- 
time receiver shown in FIG. 1 to generate the final symbol spread data signals subjected to L-path frequency-selective 
estimates b/ lc ^ for the K users. 5 fading, received at an array of P antennas comprising: 

The post-combiner interference cancellation procedure despreading with a bank of filters matched to L multipath 

described above can be continued for multiple iterations by delayed replicas of each of the K signals received at 

repeating stages 91 and 92 indefinitely. Performance gains each of the P antennas; 

due to repeated iterations diminish after the first two or three weighting the despread signals with the complex conju- 

iterations. To reduce the complexity, interference cancella- 10 gate of an estimate of the array/channel coefficient; 

tion of both types can be performed on a subset of users. summing together the weighted, despread signals to form 

What has been described is deemed to be illustrative of a K-vector; 

the principles of the invention. In particular, it should be applying a linear transform to the K-vector to suppress 

noted that the apparatus and methods of the invention may multiple access interference, said linear transform 

be implemented through various technologies, for example, ^ being the product of: 

through the use of large scale integrated circuitry, applica- the diagonal amplitude matrix (A); 

tion specific integrated circuits and/or stored program gen- the equivalent array/channel/code correlation matrix 

eral purpose or special purpose computers or (M); and 

microprocessors, including a single processor such as a the matrix inverse of the matrix (MA^'+cr 2 ]^) where 
digital signal processor (DSP), using any of a variety of is the hermitian conjugate of the array/channel/ 
computer-readable media. Further and other embodiments 20 C0( je correlation matrix, (M), is the estimated array/ 
will be apparent to those skilled in the art without, however, channel/code correlation matrix and cT is the back- 
departing from the spirit and scope of the invention. ground noise variance; and 
What is claimed is: processing the resulting signals to produce estimates of 

1. A method for detecting the data signals modulated the respective data symbols. 

respectively by multiple direct-sequence spread spectrum 25 6. The method of claim 1 wherein the linear combiner 
signals from K users, each signal being received through L employs a minimum mean-squared error criterion, 
resolvable multipath components with P receiver antennas, 7. The method of claim 6 wherein said minimum mean- 
comprising the steps of: squared error criterion designates a linear combiner for the 

a) generating matched filter outputs for at least a subset of kth user which minimizes the expected squared error 
the multiple signals, a plurality of the multipath com- 30 between the combiner output and the kth user data symbol, 
ponents of these signals, and said P plurality of receiver 8. The method of claim 1 wherein the linear combiner 
antennas; comprises a matrix W 5 formed from the product of (a) the 

b) developing coherent channel estimates from said diagonal amplitude matrix A; (b) an array/channel/combiner 
matched filter outputs; matrix which uses actual parameters rather than estimates of 

c) weighting and combining the matched filter outputs 35 the parameters; and (c) the matrix inverse of the product of 
using the coherent channel estimates; (i) the correlation matrix R; (ii) the array channel combiner 

d) employing a linear combiner which is a function of a matrix G; (iii) the square of the amplitude matrix (A 2 ) and 
correlation matrix derived from said K users 1 signals the Hermitian conjugate of the array/channel combiner 
and coefficients corresponding to said L resolvable matrix (M^) summed with the KLP-by-KLP identity matrix 
multipath components and said P receiver antennas for 40 weighted by the background noise variance o 2 , where L are 
estimating the respective data symbols for at least a the time resolvable multipath components of the signals, 
subset of the spread spectrum signals while suppressing 9. The method of claim 1 wherein the linear combiner 
multiple access interference, and comprises a matrix W fl formed from the product of (a) the 

e) passing each component of the final K-vector output of diagonal amplitude matrix A; (b) an array/channel/combiner 
the linear combiner to a decision slicer for further 45 matrix which uses estimates of the parameters; and (c) the 
processing. matrix inverse of the product of (i) the correlation matrix R; 

2. The method of claim 1 wherein said coherent channel (ii) the array channel combiner matrix G; (iii) the square of 
estimates are developed from said matched filter outputs. the amplitude matrix (A 2 ) and the Hermitian conjugate of 

3. The method of claim 1 wherein said coherent channel 5Q the array/channel combiner matrix (M^) summed with the 
estimates are developed from an auxiliary signal, KLP-by-KLP identity matrix weighted by the background 

4. A receiver comprising: noise variance cr 2 , where L are the time resolvable multipath 
a bank of correlators for despreading K code -spread data components of the signals. 

signals s^ 2 , . . . ,s^ bearing respective data b lf b 2 , . . . , 1<>- The method of claim 1 wherein the linear combiner 

b K each with at most L resolvable multipath compo- ss employs a zero-forcing error criterion, 

nents to form a K-vector; 11- The method of claim 8 wherein the linear combiner 

means for applying a linear combiner to the K-vector to forms a matrix W s from the complex block-Toeplitz corre- 

suppress multiple access interference; said linear com- lation matrix R, the complex array/channel matrix G, and the 

biner forming the product of: diagonal amplitude matrix A. 

the diagonal amplitude matrix (A); 60 12. The method of claim 8 wherein the linear combiner 

the hermitian transpose of the array/channel matrix forms a matrix W^ from the complex block-Toeplitz corre- 

(G iH ): and ^ lation matrix R, an estimate of the complex array/channel 

the matrix inverse of the matrix (RGA 2 G // +o 2 I JC£(/ ,), ma trix G, and the diagonal amplitude matrix A. 

where R is the complex block-Toeplitz correlation 13. The method of claim 1 wherein the data signals are 

matrix, cr 2 is the background noise variance, and 65 pulse amplitude modulated, and the real components of each 

^klp ^ the identity matrix whose dimension is KLPx of the combined matched filter outputs are extracted prior to 

KLP, and the linear combiner. 
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14. The method of claim 11 wherein the linear combiner 
employs a minimum mean-squared error criterion. 

15. The method of claim 12 wherein the minimum mean- 
squared error criterion designates a linear combiner for the 
kth user which minimizes the expected squared error 
between the combiner output and the kth user data symbol. 

16. The method of claim 1 wherein said linear combiner 
comprises a matrix W A formed from the product of (a) the 
diagonal amplitude matrix A; (b) equivalent array/channel/ 
code correlation matrix (M); and 

(c) the matrix inverse of the product of (i) the array/ 
channel/code correlation matrix (M); (ii) the square of 
the amplitude matrix (A 2 ) and (iii) the Hermitian 
conjugate of the array/channel/code correlation matrix 
(M ) summed with the estimated array/channel/code 
correlation matrix (&) weighted with the background 
noise variance o 2 . 

17. The method of claim 1 wherein said linear combiner 
comprises a matrix W A formed from the product of (a) the 
diagonal amplitude matrix A, (b) an estimate of the equiva- 
lent array/channel/code correlation matrix (T?I; and (c) the 
matrix inverse of the product of (i) the estimated array/ 
channel/code correlation matrix (lit), (ii) the square of the 
amplitude matrix (A 2 ) and (iii) the Hermitian conjugate of 
the array/channel/code correlation matrix (M**) summed 
with the estimated array/channel/code correlation matrix 
(M) weighted with background noise variance a 2 . 

18. The method of claim 11 wherein the linear combiner 
employs a zero -forcing error criterion. 

19. The method of claim 16 wherein the linear combiner 
zero-forcing criterion includes forming a matrix W A from 
matrix inverse of the real equivalent array/channel/code 
correlation matrix Re(M). 

20. The method of claim 16 wherein the linear combiner 
zero-forcing criterion includes forming a matrix W A from 
the matrix inverse of the real part of an estimate of the 
equivalent array/channel/code correlation matrix Re(&). 

21. A receiver comprising: 

a) space-time rake receiver having: 
i) a bank of correlators for resolving K code-spread 
data signals s a ,s 2 , . . . p K bearing respective data 
D i>t>2J • • • » b^. each with at most L resolvable 
multipath components with complex amplitude c kl 
(k»l,2, . . . ,K; 1=1,2, . . . ,L), arriving at a P element 
antenna array with complex amplitude h k jj, 
(1,2, . . . ,P) into at most KLP correlator outputs 
corresponding to at most L multipath signals for each 
of the K users at each of the P antenna array 
elements; 
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ii) means for weighting the output of each correlator 
according to the complex conjugate of an estimate of 
the coefficient for a respective channel and path 

iii) means for summing together the weighted resulting 
LP complex values for each user; and 

iv) means for forming a K-vector of the K resulting 
summations; 

b) means for applying a linear combiner to the K-vector 
to suppress multiple access interference; and 

c) means for processing the resulting signals to produce 
estimates of the respective data symbols. 

22. A receiver according to claim 21 wherein said means 
for the linear combiner comprises means for forming the 
product of: 

the diagonal amplitude matrix (A); 
the equivalent array/channel/code correlation matrix (M); 
and 

the matrix inverse of the matrix (MA 2 M /f +o 2 l5f), where 
is the hermitian conjugate of the array/channel/ 
code correlation matrix, (M) is the estimated array/ 
channel/code correlation matrix and (a 2 ) is the back- 
ground noise variance. 

23. A receiver according to claim 21 wherein the real parts 
of the K resulting summations are extracted prior to the 
linear combiner. 

24. Apparatus for use in detecting at a receiver in a 
wireless system with multiple transmitted direct-sequence 
spread spectrum signals from K users, each signal being 
received through L resolvable multipath components with P 
receiver antennas, comprising: 

a) a matched filter providing a plurality of outputs for at 
least a subset of the multiple signals, a plurality of the 
multipath components of these signals, and a plurality 
of receiver antennas; 

b) means for developing coherent channel estimates from 
the matched filter outputs; 

c) means for weighting and combining the matched filter 
outputs using the coherent channel estimates; and 

e) a linear combiner which is a function of a correlation 
matrix derived from said K users' signals and coeffi- 
cients corresponding to said L resolvable multipath 
components and said P receiver antennas for estimating 
the respective data signals for at least a subset of the 
spread spectrum signals while suppressing multiple 
access interference. 
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